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Abstract 

We investigate the effect of electric and color charge neutrality on the two-flavor color super- 
conducting (2SC) phase of cold and dense quark matter in presence of constant external magnetic 
fields and at moderate baryon densities. Within the framework of the Nambu-Jona-Lasinio (NJL) 
model, we study the inter-dependent evolution of the quark's BCS gap and constituent mass with 
increasing density and magnetic field. While confirming previous results derived for the highly 
magnetized 2SC phase with color neutrality alone, we obtain new results as a consequence of im- 
posing charge neutrality. In the charge neutral gapless 2SC phase (g2SC), a large magnetic field 
drives the color superconducting phase transition to a crossover, while the chiral phase transition is 
first order. Assuming that LOFF phases do not arise, we also obtain the Clogston-Chandrasekhar 
limit at a very large value of the magnetic field (B ~ 10 19 G) in the g2SC phase. At larger diquark- 
to-scalar coupling ratio Gd/Gs, where the 2SC phase is preferred, this limit is strongly affected by 
Shubnikov de Haas-van Alphen oscillations of the gap, indicating the transition to a domain-like 
state. 
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I. INTRODUCTION 



Color superconductivity/superfluidity (l|-|8| is now a generally accepted conjecture about 
the state of matter at high baryon densities and low temperatures. Many of its features 
have been quite comprehensively explored, and are summarized or referenced in a recent 
review The phenomenon describes Cooper pairing between quarks, with specific pair- 
ing patterns predicted as a function of "control" parameters such as external fields, quark 
masses or momentum mismatch, the latter arising usually as a consequence of neutrality 
conditions [Iol-17|. For example, at very high density fi q 3> Aqcd and for Nf = 3 massless 



flavors without any magnetic fields, the preferred pairingpattern is a flavor and color- 
democratic one termed the color-flavor-locked (CFL) phase [3j. As we tune parameters such 
as quark masses or magnetic fields to values typically found in compact stars, the CFL state 
would yield to a less symmetric pairing pattern. It is important to study these different 
manifestations of color superconductivity, since any hope of observing such exotic phases 
is tied to the physical properties of compact stars such as neutron or possibly quark stars. 
However, the analysis gets complicated due to the interplay of many effects. Such stars have 
enormous magnetic fields, of order 10 12 G at the surface, but possibly up to 10 18 G in the core, 
which affects the diquark pairing. Also, one may want to consider the partial restoration of 
chiral symmetry, since core densities are only moderately high compared to that of the CFL 
phase. Compact objects are, of course, charge neutral (either locally or globally) and color 
neutral (locally if we consider homogeneous phases) - this imposes additional constraints as 
we have mentioned. 

In this paper, we investigate the effect of electric and color charge neutrality on the 
two-flavor color superconducting (2SC) phase of cold and dense quark matter in presence 
of constant external magnetic fields and at moderate baryon densities. Several works have 
addressed parts of this problem and arrived at important conclusions: we summarize some 
of them here to orient ourselves along the new directions in this work. In the earliest works, 



based on expansions of the thermodynamic potential in powers of M s /fi, A//x [12|, |13J, it 
was found that the neutral unmagnetized 2SC phase is energetically disfavored compared 
to the CFL phase at compact star densities. However, within the NJL-model, it was de- 
termined that there are windows of the phase diagram where such a 2SC phase can arise, 
provided the diquark coupling strength (expressed as a ratio to scalar coupling) Gd/Gs is 



chosen to be large and quark masses are computed dynamically HQ. At more natural 
values of Gd/Gs, gapless quasiparticles arise when the mismatch in down and up quark 



chemical potentials 5/jl = /x e /2 > A [20|-|22j. The introduction of magnetic fields into the 
quark superconductor has important consequences. Taking the example of CFL phase, one 
of these is the emergence of a rotated photon that can penetrate superconducting quark 
matter, so that magnetic flux screening by the Meissner effect is very weak (23|. Another 
is the reduction in symmetry of the CFL phase since only d and s quarks have the same 



electromagnetic charge [24|, |25| . The magnetic field can enhance certain condensates, which 
changes the superconducting gap from the unmagnetized state [26]. Strong oscillations in 
the gap parameter and the magnetization as a function of applied field can make the ho- 



mogeneous magnetized phase unstable [27|, |28|. These works focused on the CFL phase, 



which is automatically color and charge neutral. However, neutron star densities span a 
more moderate range, where two-flavor pairing and neutrality effects play an important 



role. Electric charge neutrality of the 2SC phase with chiral condensates was studied in [29 



but in the absence of magnetic fields or the possibility of gapless phases. In previous work, 
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we studied the effect of a large magnetic field on the chiral and diquark condensates in a 
regime of moderately dense quark matter j30|. Our focus was on the inter-dependence of the 
two condensates through non-perturbative quark mass and strong coupling effects, as done 
for the unmagnetized case in [31| but neutrality effects were ignored. In such mixed 
broken phase with co-existing chiral and diquark condensate arises [3l|]. We also found that 
at large Gd/Gs, a large magnetic field changes the nature of the phase transitions in this 
mixed broken phase region. 



An extension of our work [30] to include color neutrality was performed in [32] and the 
authors found a change in the order of the phase transition from normal to superconducting 
quark matter. However, it was assumed that no intermediate phases of lower symmetry 
exist between the CFL and normal phase, which is a drastic assumption. Furthermore, 
since the NJL model has no dynamical color fields, the color chemical potential /j,s has to 
be introduced by hand, and is a very small correction to the free energy. Consequently, its 
effects on competition between phases is also likely to be small, as observed in [32|]. However, 
charge neutrality can be a more important factor for the same. In this paper, we therefore 
turn our attention to a more comprehensive analysis of the effects of neutrality (both color 
and charge) for magnetized two-flavor quark matter. This is perhaps the phenomenological 
regime most relevant to neutron stars. We present results for the quark gap, constituent mass 
and neutralizing charges as a function of density for weak and strong magnetic fields. While 
we defer the comparison of different phases to future work, this work has important new 



emergent conclusions, such as the Clogston-Chandrasekhar limit |33l. |34| for the g2SC phase, 
the disappearance of the mixed broken phase and the possibility of strongly magnetized 
domains of neutral but superconducting quark matter in the star. 

The paper is organized as follows: In section [Til we state the model NJL Lagrangian 
and its parameters. In section IHH we recast the partition function in terms of interpolating 
bosonic variables and derive the thermodynamic potential. In section IIV| we obtain the 
relevant chiral/diquark gaps as well as the neutralizing charges. In section [V] we discuss our 
numerical results for the coupled evolution of the condensates as a function of density and 
magnetic field for a fixed ratio Gp/Gs- Our concluding remarks are in section IVT1 



II. LAGRANGIAN 



We employ a NJL type Lagrangian density for two quark flavors (Nf = 2) applicable to 
scalar and pseudoscalar mesons and scalar diquarks as follows, 

£ = q H {d, ~ ieQA^ - ^T 8 G 8 ) + jlj - m] q + G s [(qqf + (01 5 fq) 2 } 

+ G D [{qi^e f e c q c ) {fi^e f 6 c q)} (1) 

where the Dirac spinor q = q ia with i = (1, 2) = (u, d) and a = (1, 2, 3) = (r, g, b) is a flavor 
doublet and a color triplet. The charge-conjugate fields of q and q are defined as q c = Cq T 
and q c = —q T C with C = — ry°7 2 . The vector r = (r 1 , r 2 , r 3 ), where the components are 
the Pauli matrices in flavor space, whereas (ef)ij and (e c ) ab3 are antisymmetric matrices in 
the flavor and color spaces, respectively. We define the chemical potential for each flavor and 
color by fi = ji — Q fi e + T 3 /i 3c + T 8 fi 8c , where Q is the generator of U(1)em, T 3 and T 8 are the 
two diagonal generators of SU (3) c , /i is the common chemical potential for non-zero baryonic 
density. Since the red and green color of a particular flavor are degenerate and the diquark 
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TABLE I: Q charges of quarks in the 2SC phase in units of e in presence of external rotated 
magnetic field B 



Flavor 


up 


down 


Color 


Red 


Green 


Blue 


Red 


Green 


Blue 


Q-charge 


+\ 


+\ 


1 


1 

2 


l 

2 






condensates in the blue color direction, we can assume /i3 C = 0. The difference of chemical 
potentials between the first two colored (red and green) quarks and the third colored (blue) 
quark of a same flavor is induced by fi$ c , and for a same color, the difference of chemical 
potentials between two flavored quarks {u and d) is induced by ji e . The electromagnetic 
charge matrix for quark defined as Q = Qf <8> l c with Qf = diag(2/3, —1/3) which coupled 
to U(l) gauge field A^. Here e is the electromagnetic charge of an electron and g is the 
SU(3) C coupling constant. The matrix rh = diag(m u , m^) is the current quark mass matrix 
in flavor basis. We take the exact isospin symmetry i.e. m u = = m . The scalar and 
diquark couplings are denoted as Gs and Gd respectively. 

We introduce four auxiliary bosonic fields which bosonize the four-fermion interactions 
via a Hubbard- Stratonovich transformation. The bosonic fields are 

a = (qq)- tx = (qi-f 5 rq) ; A = (fi^ b e f e c q) ; A* = (qil 5 e f e c q c ) (2) 

where a and n are the mesons and A and A* are the diquarks. The bosonized Lagrangian 
density becomes 

C = q [i>f (d„ - ieQAp - igT s Gl) + fa ] q-q(m + i lb Ti ■ f) q 

l AW _r \ 1 A /_ C \ cr 2 + n 2 A*A 
- g A VI H5efe c q) - - A {qil$e f e c q ) — (3) 

where m = m + a. We do not include the possibility of pion condensation for simplicity 



35[. So we set n = in our analysis. Chiral symmetry breaking and color superconductivity 
in the 2SC phase is manifest by non-vanishing vacuum expectation values (VEV) for a and 
A. In this paper we will study the effect of an external constant magnetic field B on these 
condensates. 

Since the diquark condensates of u and d quarks carry a net electromagnetic charge, there 
is a Meissner effect for ordinary magnetism, while a linear combination of photon and gluon 
leads to a rotated massless U(l) field which we identify as photon of our theory. We can 
write the Lagrangian in terms of rotated quantities using the following identity, 

eQA, + gT 8 Gl = eQA, + gf 8 G« (4) 

In the R.H.S. of the Eq. (j3J) all quantities are rotated. In flavor ® color space in units of 
the rotated charge of an electron e = y/3ge/ ^/3g 2 + e 2 the rotated charge matrix is 

Q = Q f (2) l c - 1/ (8) ^ (5) 

The generator T 3 plays no role here because the degeneracy of color 1 and 2 ensures that 
there is no long range gluon 3-field. For 2SC phase the rotated Q charges of different quarks 
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are presented in Table [B We take constant rotated background U(l) magnetic field B = Bz 
along +z axis. After integrating out the dynamical part of A M , the potential becomes 
= (0, 0, Bx, 0) in the Landau gauge. The gapped 2SC phase is Q-neutral, while overall 
charge neutrality of the matter requires a neutralizing background of strange quarks and/or 
electrons. In this article, we take the strange quark mass very large so that they do not play 
any dynamical role. 



III. PARTITION FUNCTION 

We write the partition function in the presence of a uniform magnetic field in the mean 
field approximation using the standard finite temperature path integral formalism over the 
quark fields 

Z = N J [dq] [dq}exp j jf dr J d 3 x - ^ ^ J (6) 

where N is the normalization constant, C is the rotated Lagrangian density and (3 = 1/T 
is the inverse of the temperature T. We add the kinetic energy contribution of the rotated 
U(l) external gauge field — \{F^ V ) 2 = — ^- in the partition function. The blue u and d 
quarks do not participate in the diquark condensate. Whereas, the red u paired with green 
d and green u paired with red d to form Q-charge neutral diquark condensates. Thus, we 
can write the full partition function as a product of three parts, 

Z Z C Z Ub d b Z Ur d g) Ugd r 



where Z Ub d b part is for the unpairing quarks and the Z Ur d g , Ug d r P ar t is for the quarks partic- 
ipating in pairing. The explicit form of Z Ub d b and Z Ur d g)Ug d r we present in the next section. 
The bosonized part Z c serves as a constant multiplicative factor and reads as 



Z c = Aexp <- dr d 3 x 



o 



a 2 A 2 B 2 



4G„ 4G D 2 



(8) 



The Lagrangian density C contains a diagonal 6 ® 6 chemical potential matrix fi in 
flavor eg) color space. When we incorporate the electric and color charge neutrality con- 
ditions, we have to assign different chemical potential for each flavor and color. Then the 
diagonal chemical potential matrix fi can be expressed as 

fi = diag(ii Ur , fi Ug , fi Ub , /i dr , /i dg , fi db ) (9) 

where /ij a is the chemical potential for flavor i and color a. The explicit expressions of all 
^x^ a s are as follows 

2 1 2 2 

Hu r — Vu g — I 1 — g/^e + ~^8c', AS, — \i — ~ ~H%c 

11 12 

Hd r = fid, = V + ^ + g/iSc! Vd b = V + g/ie - g/i8C (l0) 
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For the calculational simplicity, we define the mean chemical potential /x and the difference 
of the chemical potential 5 fx as 



A* = ^(AV + AHJ = t^IAS, + A*dJ = A* - g/^ + g/isc 



^ = 2^ 



AW 



2^ 



(11) 
(12) 



In order to evaluate the thermodynamic potential Q = —TlnZ/V, we introduce Nambu- 
Gorkov bispinors for each color and flavor of quark and express Z Ubdb and Z Urdg ^ Ugdr in the 
following way 



\nZ Ubdb = ihijDet^o- 1 )}; \nZ Urdg)Ugdr = ~ In (Det^ 1 )} 



(13) 



The determinant operation is carried out over the Nambu-Gorkov, color, flavor and 
momentum- frequency space and Qq 1 and Q^ 1 have the following form 



So 1 



[G 



0,Q 









6a' 



where we have used the notation A (A + ) 
read as follows 



275C/e c A(A*). The Green's functions [G^ = 



(14) 





-1 







7 M^ + eQA,±/2 7 ( 



(15) 



The determinant computation is simplified by re-expressing the Q-charges in terms 
of charge projectors in the color-flavor basis, following techniques applied for the CFL 
phase |27j]. With the color-flavor structure unraveled, we can simplify the determinant 
computation further by introducing energy projectors [3l[ and moving from position to 
momentum space using Fourier transformation, whereby we find 



InZ, 



u b d b 



{**(■ 



2/2 



E:;)} + \n{f-( P i-E-;)} 



P0,P 

+ ln 



+ In 



{/? 2 {Pi ~ E d 



(16) 



lnZ Urdg , Ugdr = 2^[ln{/3 2 (^-^:)}+ln{/3 2 (^-^!)} 



P0,P 



In ^ U-E 



(17) 



2\a\eBn 



We define E a = + p 2 + ™ 2 , if a = then pi, = pi + V \ else P ] 

for a = ±1/2, 1. Here a denotes the possible Q charges of the quarks and n labels the 
Landau levels originated in presence of magnetic field. We see from the definition of E a that 
Ei = E_i and the various other terms in the Eq. (fl6|) and Eq. ([171) are defined as 



F± 2 



E\±li u ,] E d 



Eq ± Hd b 



(18) 
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Ef = Ei± ft El = JEf + A 2 ; E± =E±±5(i 

2 2 



(19) 



The sum over po(= iu>k) denotes the discrete sum over the Matsubara frequencies. Using 
the following identity we can perform the discrete summation over p . 



£ ^ [/? 2 (Pi -E*)]=(3[E + 2Tln(l + e^ E )} ee ^ (E) 



pa 



Then go over to the 3-momentum continuum using the following replacement 

d 3 p 



(2nf 



(20) 



(21) 



where V is the thermal volume of the system. Finally, the zero-field thermodynamic potential 
can be expressed as 

* 2 A 2 „ f°° d 3 p 

\3 



m 



4G S '4G D '- J (2Tfr K)^M+^K) 
(E-) + 2 {<p (£+.) + <p (£+_) + <p (E- + ) + <p {E A _)} 



(22) 



where Q e is the contribution from non-interacting electron gas. In the limit of electron zero 
mass and small T, Q e takes the form 

■ 4 /i 2 T 2 7vr 2 T 4, 



12tt 2 



6 



+ 



180 



(23) 



In presence of a quantizing magnetic field, discrete Landau levels suggest the following 
replacement 



d 3 ~p \a\e 



(2tt)' 



<,7T Z 



POO 

^ OL n \ dp z where a n = 2 - <5 n0 

n=0 



(24) 



where a n is the degeneracy factor of the Landau levels (all levels are doubly degenerate 
except the zeroth Landau level). The thermodynamic potential in presence of a magnetic 
field is given by 

,2 A 2 



= — — I" TFi 1" + 5^ where a 6 0,1,- 



4G 5 4G D 

The Q a is the contribution to the Qb from the quarks of Q-charge a and they are 

d 3 p 



(25) 



n 



o (2tt) c 



(26) 



eB 



n=0 



(27) 



Hi 

2 



57T 



2 Z^"" 

n=0 



dp z 



(E+ + )+0(E+_)+0(E A+ )+0(E A _) 



(28) 
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IV. GAP EQUATIONS AND NEUTRALITY CONDITIONS 



Minimizing the thermodynamic potential in Eq. (J25]) w.r.t. m and A, we can derive 
chiral and diquark gap equations respectively. We can also get the color and electric charge 
neutrality conditions minimizing the thermodynamic potential w.r.t. fi 8c and /i e respectively. 



A. Chiral gap equation 

The gap equation for chiral condensate can be obtained from the following equation 

dtt B (a, A, /i 8c , /i e ; /i, eB, T) 



da 







(29) 



dQ 



B 



a 



^^{i -/(*!) -/(*£)} 



da AGs '[Jo (27r) 3 £ 

~D 00 POO 1 

„_n J —oo 



n=0 
oo 



pD ^ POO i p- 



16vr 2 



n=0 



where /(x) = [1 + exp 1 is the Fermi-Dirac distribution function. 



(30) 



B. Diquark gap equation 

The gap equation for diquark condensate can be obtained from the following equation 

dQ B (a, A, /i 8c , fj, e ; y, eB, T) 



dA 







(31) 



dQ 



B 



dA 



A 



1 



IB 



AG D 16tt : 



n=0 



^2® n / dp 



2^{1 -f(E- + )-f(E-_)} 



+ 2-^{l -/(£++)- /(£+_)} 

-C/A 



(32) 



C. Color charge neutrality 

The total system should be color neutral. We can impose the color charge neutrality 
condition by choosing /x 8c in such a way that the system has net color charge n 8c zero. It 
follows that, 



dQ B (er, A, use, He] H, eB, T) 



dfx 







(33) 



8c 



s 



9/i 



8c 



dp 2 



(2tt 
16^2 



n=0 



-/(£"+)- /(£"_)} 

A 



(34) 



D. Electric charge neutrality 

The total system should be electrically neutral. We can impose the electric charge neu- 
trality condition by choosing \i e in such a way that the system has net electric charge n e 
zero. It follows that, 



dVt B (a, A, fi 8c , fi, eB, T) 







(35) 



on 



B 



d/j, e 



+ 



o 

eB 

167T 



_^[2{/(Ej)-/(^)}] + ^£«n / ^W/(K)-/(^)}] 



n=0 



£ 

n=0 



a, 



dp 2 



2^{1 - f(E A+ ) - f(E-_)} - 2^{1 - f(E+ + ] 
E A & A 



- /(£+_)} - 6{f(E A+ ) + /(£++) - /(£+_) - f(E A .)} 



+ ^§ = 



(36) 



We derive two gap equations and two charge neutrality conditions above for finite T, but 
our all numerical results in the following section is in the T — > limit. We investigate the 
finite T results elsewhere. 



V. NUMERICAL ANALYSIS 



In order to assess the effects of electric and color charge neutrality on quark matter in 
presence of an external magnetic field, we solve the two gap equations with electric and 
color charge neutrality conditions imposed. Since this set of four equations all involve inte- 
grals that diverge in the ultra-violet, we regularize in order to obtain physically meaningful 
results. It is common practice to use a sharp cutoff 1 to regulate, but this can introduce 
abnormal behavior in many thermodynamical quantities, especially when a system with a 
discrete energy spectrum is considered. In this numerical study, we therefore use a smooth 
cutoff function while performing the numerical integrations in momentum space. There is a 
variety of smooth regulators used in the literature |3|, |27|, and all are expected to yield 
qualitatively similar results. We use a Fermi-Dirac type of smooth regulator (26| 



fc{Pa) 



tanh 



Pa 



A 



a 



(37) 



1 An exam 



iple is the step function 9 (A — |p|) in momentum space where A is the cutoff momentum. 



9 



where p a = \Jp\ a + V% if a — then p 2 ±0 = p x +p 2 else p 2 ± a = 2\a\eBn for a = 1, ±1/2. In 

Eq. (1571) A is the momentum cutoff and a is a free parameter that measures the sharpness of 
the cutoff scheme. Larger values of a imply that the cutoff is gentler. We choose a = 0.01A 
throughout our numerical analysis. We fix a momentum cutoff A and NJL model parameters 
in the chiral limit tuq = as given in j29| . 



A = 0.6533 GeV, G s = 5.0163 GeV" 2 and G D = pG s (38) 

where p is a free parameter. Fierz transforming one gluon exchange implies p = 0.75 for 
number of colors N c = 3 and fitting the vacuum baryon mass gives p = 2.26/3 [37]] . Since the 
underlying interaction at moderate density is bound to be more complicated, we therefore 
choose to vary the strength of the diquark channel to investigate the effect of p on the 
magnetic color superconducting phase. In [2l[ where no magnetic field is considered, it is 
found that the ground state of charge neutral two flavor quark matter is very sensitive to 
Gd- In particular, one finds the 2SC phase (i.e. A > 5p) for p > 0.8, the g2SC phase (i.e. 
A < 5p) for 0.7 < p < 0.8 and the normal quark matter (i.e. A = 0) for p < 0.7. We take 
the bare quark mass m = throughout our analysis. One can consider realistic nonzero 
m and get a set of model parameters by fitting three vacuum quantities, namely the pion 
mass m^, the pion decay constant /„- and the constituent quark mass m(p = 0). For nonzero 
magnetic field, the model parameters will differ only slightly from zero field values. For our 
analysis, we keep model parameters fixed at zero field values. We note that for the choice 
of a very sharp cutoff function, Eq. (124"]) becomes 

-A' r a 2 i 

(39) 



f°° d 3 p , , , f A d 3 p \a\eB n ^ f A ' T A 

/ T^^P) ~ / 7T^3 ~3Zr 1^ a n / d Pz where n max = Int 

^0 (27T) Jo (2-R) 07T ^ J_ A , 



n=0 



2aeB 



and A' = a/A 2 — 2\a\eBn. Using the fact that p 2 z > 0, we can find the maximum number of 
completely occupied Landau levels n max . In the limit of weak magnetic field n max becomes 
very large and the summation over discrete landau levels becomes quasi-continuous. The 
zero magnetic field results remain an excellent approximation in the weak magnetic field 
limit, as we see below. 



We reiterate that our main aim in this work is to investigate the effect of magnetic fields 
on the chiral (m) and diquark (A) condensates for different pairing strengths (p) and with 
the electric and the color charge neutrality conditions imposed. Whether the quark matter 
is in g2SC phase or 2SC phase depends on the value of p. Therefore, in this work we take 
three different values of p=0.75, 1.15 and 1.4. We may consider the effect of electric and 
color charge neutrality on the chiral and diquark condensates for the following four separate 
cases, of which the last two are new : (I) no charge neutrality condition, which was already 
discussed inj30j; (II) only the color charge neutrality condition, already discussed in[32{]; (III) 
only the electric charge neutrality condition is imposed; (IV) both the neutrality conditions 
are added. 

Fig. 0Ja) shows m, A and the electron chemical potential p e in neutral two-flavor quark 
matter for p = 0.75 in the weak magnetic field limit (eB = 0.005 GeV 2 ) - this is the g2SC 
phase. Fig. []Jb) shows the small color chemical potential /i§ c . We use a sharp cutoff for 
the results in Fig. [1] in order to compare with the zero-field results discussed in [29| . For 
a magnetic field eB = 0.005 GeV 2 , n max is of the order of 50 and the discrete summation 
becomes almost continuous. We find that our results for m, A, p e and p$c in the case of 
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FIG. 1: For electric and color charge neutral quark matter with small magnetic field: (a) The 
dependence of m, A and p e on p for eB = 0.005 GeV 2 taking p = 0.75 and (b) the corresponding 
dependence of ps c on p. Here we use a step function cutoff. 



eB = 0.005 GeV 2 are in very good agreement with the zero magnetic field results. In Fig. 
QJa) we see that for p, < p, c (« 0.315 GeV) the quark possesses a large constituent mass 
m (chiral symmetry broken) and in the \x > \x c region, chiral symmetry is restored. This 
phenomenon is first order in nature. With p, c as the critical chemical potential for this phase 
transition, the p, < p, c region is the chiral symmetry broken (CSB) phase. In the p, > p c 
region, the diquark condensate appears, and this is the color superconducting (CSC) phase. 
We find here that p c ~ m(fi < fi c ) and we have checked that this is true for other sets 
of parameters as well. In the rest of our analysis wherever it is not mentioned we use the 
smooth cutoff from Eq. [37] with a = 0.01 A. 

In Fig. [2] we show several physical quantities in neutral quark matter for large magnetic 
field eB = 0.1 GeV 2 , for two different values of p. Figs. [2(a) and [2(b) show the behavior of 
m, A, 8p, and A as a function of p, for p = 0.75 and p = 1.15 respectively. Here, A is the 
diquark gap when no charge neutrality condition is considered. In Fig. [2(a) 5p > A and 
thus, this is the g2SC phase. The value of A in the neutral g2SC phase is much smaller than 
the A . Observe that in Fig. [1(a), the chiral and CSC phase transition occur at the same 
point p c and both transitions are first order in nature. Whereas in Figj2j(a), there is a small 
window between the chiral and CSC phase transition. Though the chiral phase transition 
is still first order in nature, CSC phase transition to the g2SC phase becomes a crossover in 
presence of magnetic field. We see in Fig. [2(b) that 5p> < A and thus, this is the 2SC phase. 
In contrast to Fig. [2(a), A in the neutral 2SC phase is not that much different from A . 

In Fig. [3]we show the effects of choosing different cutoff schemes on m and A. In Fig. 
[3(a), we show the effect of cutoff on m as a function of eB in the CSB phase at p, — 0.2 
GeV taking p = 0.75 with no neutrality condition is imposed (case I). In Fig. [3(b), we show 
a similar plot for A in the CSC region taking p, = 0.4 GeV. We choose /i = 0.2 GeV from 
the region p, < /i c (CSB phase) and /i = 0.4 GeV from the region p, > p, c (CSC phase) 
as the representative quark chemical potential in these phases. For these plots we use the 
smooth cutoff function shown in Eq. ([37j) for three different values of sharpness parameter 
a and also a sharp cutoff. In all cases, we observe an oscillatory behavior in m and A as 
a function of eB. The ripples in the gap parameters become steeper for smaller values of 
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(a) p = 0.75, eB = 0.1 GeV 2 



(b) p = 1.15, §B = 0.1 GeV 2 
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FIG. 2: For electric and color charge neutral quark matter with large magnetic field: (a) The 



dependence of m, A, Sfi and A on p for eB 
p = 1.15 (right panel). 



0.1 GeV 2 taking p = 0.75 (left panel) and (b) 



(a) p = 0.75, \i = 0.2 GeV (b) p = 0.75, \i = 0.4 GeV 
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FIG. 3: The effect of sharp/smooth cutoff on the chiral and the diquark gap: (a) The dependence 
of m on eB at p = 0.2 GeV in the chirally broken phase and (b) The dependence of A on eB at 
p = 0.4 GeV in the superconducting phase. In both cases, neutrality conditions are not imposed. 
A smooth cutoff given in Eq. (I37p with a = 0.1 A, 0.05A, 0.01 A and a step function cutoff are used. 



a. For the step function the ripples reach its maximum height. We find a dip in the curves 
when A 2 / (2\a\e~B) takes an integer value. This kind of oscillatory behavior is analogous to 
the well known Shubnikov de Haas-van Alphen oscillation observed in metals in presence of 
magnetic field at very low temperature. Due to the Landau level quantization in presence of 
magnetic field the density of states and hence, various thermodynamic quantities oscillate 
as a function of magnetic field. The wavy behaviour of m and A diminishes as we go to 
smaller values of eB. Finally m and A tend toward their zero field values for sufficiently 
small eB which corresponds to the quasi-continuous limit. We checked for different values of 
p and obtained similar oscillatory behavior for the gaps. Since only Gs enters in the chiral 
gap equation, m remains same for different values of p. On the other hand, A increases (but 
keeps the same oscillatory features) as we increase p (or equivalently the diquark coupling 
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FIG. 4: (a) The dependence of A on eB at fi = 0.4 GeV in the CSC phase for case I and case II 
(see text) taking p = 0.75 and (b) the dependence of corresponding /ig c on eB for the case II. 
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FIG. 5: (a) (top and bottom panels) The dependence of A and p, e on e.B at p = 0.4 GeV in the 
CSC phase for the case III and case IV (see text) taking p = 0.75 and (b) the dependence of 
corresponding p$ c on eB for case IV. 



strength Go) since Ge> is directly involved in the diquark gap equation. 

In Fig. 11(a) we show A as a function eB in the CSC phase for case I and case II at 
/i = 0.4 GeV taking p = 0.75. We see that the effect of color neutrality (/isc) on the 
CSC phase is very small. Thus, we can sometimes neglect fi$ c for computational simplicity 
without changing the conclusions qualitatively. In Fig. H(b), we show /ig c as a function of 
eB corresponding to case II of Fig. H(a). In j2£j it was shown that fi$ c can be positive or 
negative as a function of /i. Here we find similar behavior in /i 8c as a function of eB for a 
fixed value of /i = 0.4 GeV. In the weak field limit, the value of /ig c is only about a few MeV, 
but its value can be fairly large for large values of eB. Physically, the magnetic field creats 
a stress between the Fermi surfaces of the blue and the red/green quarks. 

In Fig. |5(a) we show the behavior of A and /i e in the g2SC phase (p=0.75) as a function 
of eB taking pi = 0.4 GeV. We compare the results without and with color charge neutrality 
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(a) p = 0.75, \i = 0.4 GeV ( b ) p = .75, \i = 0.4 GeV 




FIG. 6: (a) The dependence of 6(i° (= pP e /2) and A°/V2 on eB at fi = 0.4 GeV (b) (top panel) a 
comparison between /i e (= /i®) in the normal quark matter and \i e in the CSC phase as a function 
of eB at fi = 0.4 GeV without ^sc (b) (bottom panel) a comparison between A, 5[i and A°/\/2 as 
a function of ei? at \x = 0.4 GeV without ^§ c . We take p = 0.75 for both the plots. 




FIG. 7: (a) The dependence of A, 5fi, S/jP and A°/V2 on eB at fi = 0.4 GeV for case IV (see 
text) taking p = 1.15. (b) The dependence of A and 5ji on eB at /U = 0.4 GeV for case IV taking 
p= 1.4. 



condition imposed. Fig. [5](b) shows the corresponding /j, 8c for the bottom panel of Fig. [5]^a). 
The effect of /i 8c on A and /i e is small, as before. But if we compare the diquark condensate 
A in Fig. Ufa) with that in Fig. IHa), we see that A in the charged CSC phase has a 
much larger value than the neutral case. This is because in the g2SC phase the presence of 
/i e creates difference in the Fermi surfaces of the two pairing quarks and reduces the value 
of the diquark gap. Also in contrast to Fig. H(a), enforcing charge neutrality drives A to 
zero at a certain value of eB. It was argued in [Io[ that BCS pairing is not possible if 
S/jP > A°/^/2. This is clearly confirmed from Fig. [6]^a) for the g2SC phase, where we see 
that Sfi°(= /ig/2 where is the electron chemical potential in the normal quark matter) is 
always less than A / \/2 until A (the gap in the neutral phase) disappears, at which point 
SfjP beats A / \/2. From Fig. Mjo) (top panel), we also see that [i e is larger in the CSC phase 
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than that in neutral normal quark matter. At around eB ~ 0.11 GeV 2 , A becomes zero 
and p e in the CSC phase coincides with the p e in the neutral normal quark matter. Thus, 
the magnetic field is playing an important role in delineating the criterion for breakdown 
of the homogeneously paired phase. At this critical magnetic field, an alternate phase, 
such as the LOFF phase would be preferred. This is the familiar Clogston-Chandrasekhar 
limit predicted for electronic superconductors. An important difference is that in condensed 
matter systems, there is no analog of the charge neutrality condition, so A would go to 
zero. In case of quark matter, where QCD provides the pairing force, it is the charge neutral 
gap and not A that vanishes. 

Finally, we show in Fig. [7] the behavior of A and Sp at p = 0.4 GeV for two large values 
of p = 1.15 (left panel) and 1.4 (right panel). In presence of a magnetic field, we observe 
that even if p > 0.8, quark matter is not necessarily in the 2SC phase. In Fig. [3(a) we 
see that most of the time A is larger than Sp which is the 2SC phase but in the range of 
eB ~ 0.4 — 0.5 GeV 2 , A becomes zero. This happens again at the Clogston-Chandrasekhar 
limit, but at a magnetic field value that is about 4 times larger than for the g2SC phase. 
This is because the mismatch Sp in the charged CSC phase does not get as severely stressed 
by the magnetic field. We also note that due to the vanishing gap, this phase can no longer 
be classified as the 2SC phase. In Fig. [7(b) we see that for very large values p = 1.4, A is 
always above Sp in the entire range of eB. Thus, this phase is the usual 2SC phase. 



VI. CONCLUSIONS 

We have studied the effect of a large magnetic field and neutrality constraints on the chiral 
and diquark condensate in a two-flavor superconductor using the N JL model. We used a self- 
consistent scheme to determine the condensates, numerically iterating the coupled (integral) 
equations for the chiral and superconducting gap under the constraints of color and charge 
neutrality. We presented results for the 2SC as well as the gapless 2SC phase by choosing 
various values of the Gd/Gs ratio. For magnetic fields B < 10 18 G, our results for the chiral 
and superconducting gaps align with the zero magnetic field results of j29|, since many 
Landau levels (~ 50) are being summed over, which approximates the continuum result. 
We note that the chiral and superconducting phase transitions occur at the same value 
of p q . This changes at very large values of B > 10 18 G for the gapless 2SC phase, where 
the color superconducting phase transition becomes a crossover. Furthermore, the value of 
the superconducting gap in the gapless 2SC phase is much smaller than that in the charge 
neutral 2SC phase. Both effects have to do with stresses on the up and down quark Fermi 
surfaces due to the charge neutrality constraint. 

As found in previous works on the three-flavor gap or CFL phase 27|, 28[, the two-flavor 



gap at large magnetic field B > 10 18 G exhibits Shubnikov de Haas-van Alphen oscillations 
due to discrete fluctuations of the quark density of states as a function of magnetic field. 
These oscillations can become quite large for strong magnetic fields, and even for Gd/Gs 
ratios that normally correspond to strong coupling, i.e, yield the 2SC phase at zero field, 
the gap can vanish when 5p > A. Physically, this means that the homogeneous gap ansatz 



is not the correct one anymore, and probably a LOFF phase |38|, |39| is the preferred one. 
Large oscillations in the gap can also lead to the formation of magnetic domains if the ho- 
mogeneous phase becomes unstable/metastable, with applications to magnetar phenomena 



as argued in [27|. Another interesting consequence of charge neutrality, through the mo- 
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mentum mismatch of the down and up quarks, is the Clogston-Chandrasekhar limit |33l. l34j. 
originally formulated as a critical magnetic field in hard electronic superconductors. While 
difficult to confirm experimentally in metals, due to the strong Meissner effect, it can be 
achieved more easily in a color superconductor, where the rotated photon 7 can penetrate. 
In the gapless 2SC phase, we find that this limit is reached at B pa 2 x 10 19 G. It should 
be noted, however, that it is not the magnetic field that destroys the BCS state. In fact, 
without electric charge neutrality enforced, the magnetic field would actually increase the 
value of the gap, since the magnetic moments of the pairing quarks are aligned. Rather, the 
mismatch of Fermi surfaces increases with increasing magnetic field, forcing the gap in the 
(locally) charge neutral phase to get smaller. Thus, the magnetic field has an essential role 
to play in delineating the conventional Clogston-Chandrasekhar limit 5/jP = A / \/2. 

Finally, although we have used an effective model with a cutoff, which causes some 
uncertainty in the calculated values of the quark gap and constituent mass, it is interesting 
to speculate on possible physical consequences of the large magnetic field, and its effect on 
the condensates. At some fixed large value of the local field B and in the density window 
of the metastable regions mentioned above, cooling of two-flavor quark matter below the 
superconducting transition temperature T c can result in the formation of domains or nuggets 
of superconducting regions with different values for the gap, even at constant pressure (a 
detailed analysis requires assessment of surface and screening effects in superconducting 



quark matter (40 , HI). We note that once neutrality is imposed, there is no range of density 



where the chiral and diquark gap co-exist. The mixed broken phase found in [3JJ does not 
survive under physical conditions of neutrality. However, in the chirally broken phase (where 
the diquark gap is zero), due to the magnetic field, even the chiral gap (or constituent mass) 
oscillates. We can also imagine that external magnetic fields of order eS//i 2 ~ 1 show 
some local variation on the microscopic scale in the initial stages of the formation of the 
dense neutron star. In this case, magnetic domains with different magnetization can form. 
At lower density, domains of broken chiral symmetry can also form, hearkening back to 



the disordered chiral condensate (DCC) idea [42|. Such kinds of nucleation and domain 
formation will release latent heat that might be very large owing to the large value of the 
magnetic field, serving as a large internal engine for possible energetic events on the surface 
of the neutron star. Such internal mechanisms are unlikely to occur in a pure neutron star 
without a quark core 43j and could have important applications for astrophysics of gamma- 
ray bursts [44| . In addition, it has been recently shown that large magnetic fields in strange 
quark matter introduce anisotropies in bulk viscosities, which can change the stability region 



of the r-mode in such stars [45| and the epoch of gravitational wave emission as the star spins 
down. A similar effect may arise in magnetized superconducting two-flavor quark matter, 
pointing to another important link between the microscopic and astrophysical consequences 
of a large magnetic field in quark matter in the core of neutron stars. 
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